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Abstract 

Physics of dense matter is an extremely rich and complex scientific field, 
mixing laboratory experiments and astronomical observations with complex 
calculations. The aim of this lecture is to discuss in some detail two simple 
but important problems which have considerable importance for the devel- 
opement of theoretical work in semiclassical dense matter physics. 

Introduction 

The first historically known attempt to study experimentally a material 
under high presssure dates from the XVIII century. An english "gentle- 
man" named Mr. Canton tried to compress water to a pressure of the order of 
0.1 GPa (Block et.al.,1980).To his astonishement,he observed a "transforma- 
tion" of water into ice. In the vocabulary of modern condensed matter physics 
this would be called a phase transition. For a recent example of a study of 
water and ice under high pressure see (Sundberg and Lazor,2004).We can 
only wonder what would Mr. Canton say todaywhen 15 different phases of 
water ice are known. Systematic laboratory high pressure studes started 
near the end of the XIX century owing to the efforts of P.W.Bridgman at 
Harvard University (Bridgman, 1964). High pressure experiments have been 
revolutionized by the discovery of the diamond anvil cell (for a recent review 
see for example Struzhkin,Hemley and Mao, 2004). Experimental and theo- 
retical research in high pressure (or, more generally, dense matter physics) is 
pursued throughout the world. 

Scientific research in high pressure physics has started in Serbia at the 
beginning of the sixties (Savic,1961).In that paper, which was short and sim- 
ple but initiated an important series of publications, he pointed out an in- 
teresting fact: namely,the mean densities of the planets, as derived from 
observational data known at the time, could be related to the mean solar 
density by an extremely simple equation. He did not make any hypothesis 
about the physical nature of this correlation. 



In the next four years Savic in collaboration with Kasanin (Savic and 
Kasanin, 1962/65) developed a theory of the behaviour of materials under 
high pressure.lt was later nick-named the SK theory for short;another syn- 
onim is "the semiclassical theory" . This last term is justified by the fact that 
this theory uses only the basic notions of atomic and molecular physics,while 
it is largely founded on laws of classical physics. 

The value of this theory is that it offers the possibility of preforming in 
a simple way calculations which would be prohibitively complex in rigor- 
ous statistical mechanics. Logically,this simplicity very often diminishes the 
precision of the results. 

The basic physical idea of their theory is that subduing a material to 
sufficiently high pressure leads ultimately to excitation and ionisation of the 
atoms and molecules which constitute it. 

At the time when they were starting to develop their theorythe study of 
the influence of external fields on atoms and molecules was not a theoretically 
well developed subject. However, Savic and Kasanin used experimental data 
known at the time to support their idea. 

NamelyBridgmann's experiments were showing that in some materials 
under certain values of external pressure first order phase transitions oc- 
cur. It was also known from the analysis of the propagation of seismic waves 
that their speed changes abruptly a certain depths in the interior of the 
Earth. Details on their theory are readily available in the literature (for 
example, Celebonovic, 1995) . 

In the remainder of this contribution,we shall discuss at some length two 
problems of considerable importance for the foundations of the SK theory: 

the relationship between the mean planetary densities derived from ob- 
servation and the mean solar densityand 

the calculation of the phase transition points in this theory. 

The density problem 

The basic aim of (Savic, 1961) was to derive some form of relationship be- 
tween the solar density and the mean planetary densities. The densities were 
calculated from the values of planetary masses and radii known at the time. It 
was shown in that paper that it was possible to fit the values of p by the 
following simple expression: 

P = Po2 4 ' (1) 

where po = 4/3 is the value of the mean solar density known at the time 
expressed in CGS units. By suitably assigning to the exponent (j) various 
integral values it turned out to be possible to reproduce known values 
of the planetary densities. This fit is extremely important for the complete 
developement of the SK theory because it determines the ratio of densities 
in two successive phases of a material under high pressure. However, as 
eq.(l) was derived more than 40 years ago,the logical question about the 
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Table 1: Densities of planets and satellites 



object 


Pmean 





Pcalc 


Sun 


1408 





1408 


Mercury 


5427 


2 


5632 


Venus 


5243 


2 


5632 


Earth 


5515 


2 


5632 


Moon 


3350 


5/4 


3349 


Mars 


3933 


3/2 


3982 


Phobos 


1900 


2/5 


1858 


Deimos 


1750 


1/3 


1774 


Jupiter 


1326 


-1/10 


1314 


JI 


3530 


7/5 


3716 


JII 


3010 


1 


2816 


JIII 


1940 


1/2 


1991 


JIV 


1840 


2/5 


1858 


Saturn 


687 


-1 


704 


Titan 


1881 


2/5 


1858 


Uranus 


1270 


-1/7 


1275 


Ariel 


1670 


1/4 


1674 


Neptune 


1638 


1/4 


1674 


Triton 


2050 


1/2 


1991 


Pluto 


1750 


1/3 


1774 


Charon 


2000 


1/2 


1991 



possibilities of its improvement or even change has started to appear. In the 
meantime the values of planetary and solar densities have been improved, so 
this also presented a challenge to eq.(l). 

The following table contains modern values of the mean densities of 
various objects within the Solar System. All the values given in the table 
were taken from the internet address http:/ /nssdc. gsfc.nasa.gov/planetary/ 
planetary home.html . The column denoted by p mea n contains densities 
of all the planets of the Solar System and their most important satellites 
derived from their observed masses and diameters. The values of the 4> in the 
third column of table 1 have been obtained by rounding off values calculated 
from eq.(l). 

The column with the heading p ca i c contains values of the density calcu- 
lated by eq.(l) and the values of <j) given in the table. 

Several conclusions can be drawn by simple inspection of Table 1. 

The aim of this work has been to test eq.(l) on modern data and on 
a larger sample of bodies than the one used in (Savic,1961). Eq.(l) pro- 
posed by Savic and Kasanin more than 40 years ago has been confirmed 
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on a much larger sample of objects. Note that non-integral values of the 
exponent <j) occur in Table 1. This marks a considerable difference with 
the original paper by Savic, in which only integral values of <p occured.A 
possible interpretation of this result can perhaps be found in considerations 
of the similarities between atomic sturcture and the structure of the plan- 
etary system. The physical meaning of eq.(l) is still a problem open for 
discussion, but the important concusion of the present work is that the fit 
proposed by Savic and Kasanin is valid for a much larger set of objects, under 
the condition that the exponent 4> takes on non-integral values. 

The phase transition points 

The theory proposed by Savic and Kasanin predicts that a material 
subdued to increasing pressure undergoes a sequence of first order phase 
transitions, numbered by an integer index i.The algorithm for the calcula- 
tion of these values of pressure is presented in detail in previous papers on 
the subject (for example Celebonovic,1995 or later publications). Just for 
illustration,the final expression for the phase transition pressure is 

p tr = 0.5101p*;z = 1,3,5,.. (2) 

and 

Ptr = 0.6785p*;z = 2,4,6... (3) 

where 

p* 1.8077 f3i(V)- 4/3 2 4i / 3 M Bar (4) 

where V is the molar volume of the material under study,and the symbol 
Pi denotes a constant determined within the theory for each phase of the 
material. 

However, the problem with eq.(2),and in general with various calculations 
within this theory is that they have a weak contact with the well established 
theories of statistical physics. In the remaining part of this paper one such 
connection will be proposed. 

The starting point of all rigorous considerations in statistical physics is 
the Hamiltonian of the system under consideration. Starting from it, one can 
derive the expression for the partition function, and from it all the thermo- 
dynamic potentials. Accordingly,in attempting to link the theory proposed 
by Savic and Kasanin with modern statistical physics, one should start from 
a Hamiltonian (or a suitable expression for one of the thermodynamic po- 
tentials) which can in principle predict a number of phase transitions of 
first order. The order of the transitions is important here because the the- 
ory we are discussing here is limited to first order transitions. In the Landau 
theory of phase transitions,the expression for the Gibbs potential per unit 
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volume of a system having spontantenous magnetization M has the form 
(for example,LeBellac,1988): 



g(M) = ^r (T)M 2 + ±u M 4 + ±v M 6 (5) 

where uq < 0, vq > and ro has the form r$(T) = a(T — To) ; To denotes 
some critical temperature relevant for the problem under consideration and 
a is a constant. 

It can be shown that g(M) = for 

M 2 = _ 15no [900ng - lUOrpVo} 1 / 2 

v 2v U 

Clearly,this expression is real if no < 0,vq > and uq > [(f)^o^o] 1 ^ 2 - 
It follows from eq.(4) that dg/dM = for 

M 2_ IQnp t [400^-480r z;o] 1 / 2 
v 2v 

which implies that Uq > (|)rofo- In a similar way one could determine 
the position of the points in which the second derivative of g{M) becomes 
equal to zero. 

Replacing the magnetization M in eqs.(3)-(5) by the pressure P renders 
these expressions applicable to the case of materials under high pressure. In 
this way a link between the theory proposed by Savic and Kasanin and 
the Landau theory of phase transitions becomes possible. Eq. (2) contains an 
experimentally known parameter (V) and another one derived within the 
theory (/?). 

Combining eqs.(2),(3) and (6) and multiplying out all the pure num- 
bers, it follows that for i = 1: 

_ 15no _ [900^^-1440^ = 2(F) - 8/3 (g) 

vo 2v ' iy 1 y 1 

where value of the constant f3\ is known from the theory of Savic and 
Kasanin. This means that eq.(8) gives a direct connection between this theory 
and the Landau theory of phase transitions. Reasoning in the same way for 
i = 2,one could obtain a similar equation. Solving a system of these two 
equations one could obtain the values of the parameters of the Landau theory 
expressed in terms of the parameters of the theory of Savic and Kasanin. 
This calculation is algebraically tedious but straightforward. As a result one 
would obtain expressions for uo and vo in terms of ro(T) and the parameters 
which enter the calculations within the SK theory. 
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Conclusions 



In this lecture two apparently simple problems in dense matter physics 
have been addressed to some extent. Both of the problems discussed are 
related to the foundations of the theory of behaviour of materials under 
high pressure proposed by P.Savic and R.Kasanin nearly fifty years ago. 
The general conclusion is that one of the basic statements of the SK theory 
is validated with modern data, and that this theory can be connected with 
modern theories of phase transitions. 

The simple expression proposed in (Savic,1961) which links the mean 
planetary densities with the mean solar density was tested on modern data. 
All the planets and several satellites of each of them were included -21 
object in all. It was shown that the original fit (eq.(l) of the present paper) 
is confirmed with modern data. However, the exponent (ft in eq.(l) is no longer 
integral but can also take non integral values. This conclusion, although 
important for the SK theoryopens up the problem of the interpretation of 
non-integral values of the exponent (ft. 

In order to make a contact between the SK theory and some of the 
better known theories of phase transitions, it was combined with the Landau 
theory of first order phase transitions. The result is that the Landau theory 
also has several first order phase transition points, and that the parameters 
of the Landau theory can be expressed in terms of the parameters of the SK 
theory. 
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